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ABSTRACT 


This  report  is  the  first  volume  of  the  RADC  Method  of  Moments  Appli- 
cations Series.  The  volumes  of  this  series  deal  with  the  application  of 
the  method  of  moments  [l]-[2]  to  various  practical  problems  such  as  antenna 
coupling,  near  field  prediction,  aperture  coupling,  feedline  interference, 
etc.  These  problems  are  of  particular  interest  to  antenna  system  designers 
and  those  in  the  area  of  electromagnetic  compatibility  [3]-[4].  There  is 
an  emphasis  in  the  Applications  Series  upon  experimental  verification  of 

effects  noted  and  cross  checking  of  results. 

\ 

^hie  first  volume  serves  the  purpose  of  an  introduction  to  the  Method 
of  Moments  Applications  Series.  The  basis  of  the  method  of  moments  is  out- 
lined and  the  notation  and  conventions  used  throughout  the  series  are  intro- 
duced. Illustrative  examples  are  presented. 

The  following  volumes  of  the  series  have  been  completed  at  this  time; 
several  other  volumes  are  in  preparation: 

VOLUME  1 - AN  INTRODUCTION  TO  THE  METHOD  OF  MOMENTS 

VOLUME  2 - NEAR  FIELDS  OF  THIN-WIRE  ANTENNAS 

VOLUME  3 - THIN  WIRE  ANALYSIS  PROGRAM  (TWAP)  USER  MANUAL 


VOLUME  4 - FEEDLINE  INTERFERENCE  WITH  DIPOLE  PERFORMANCE  - THEORY 
AND  EXPERIMENT 


VOLUME  5 - THIN-WIRE  ANTENNA  ANALYSIS  COMPUTER  CODES  COMPARED 
WITH  MEASURED  DATA 


VOLUME  6 - MATRIX  METHODS  FOR  STATIC  MICROSTRIP  PROBLEMS 
VOLUME  7 - INTERNAL  SYRACUSE  UNIVERSITY  MEMORANDA 
VOLUME  8 - APERTURE  COUPLING  THROUGH  LONG  SLOTS 


VOLUME  9 - NEAR  AND  FAR  FIELD  THIN-WIRE  COUPLING 
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7.  - THE  BASIS  OF  THE  METHOD  OF  MOMENTS 


1.1  INTRODUCTION 

There  have  been  several  dramatic  shifts  in  the  direction  of  technical 
fields  over  the  last  few  decades.  This  is  especially  true  of  the  field  of 
electromagnetic  theorv  and  its  application  to  radiation  and  scattering  pro- 
blems. During  the  last  few  decades,  this  field  has  evolved  continuously 
from  one  based  solidly  on  the  classical  theory,  but  applicable  to  a narrow 
range  of  problems,  to  one  based  on  the  classical  theory,  but  formulated  in 
terms  suitable  to  the  computer  age  and  applicable  to  a wide  variety  of 
practical  problems. 

Electromagnetic  problems  can  be  formulated  in  two  alternative  ways 
either  as  (1)  a partial  differential  equation  with  boundary  conditions, or 
as  (2)  an  integral  equation.  During  the  1930's  and  1940's,  emphasis  was 
placed  upon  the  corresponding  classical  methods  of  separation  of  variables 
and  the  classical  integral  equation  solutions.  A high  degree  of  ingenuity 
was  used  to  apply  these  classical  methods  to  practical  problems.  However, 
only  a narrow  range  of  problems  could  be  tr sated  bv  these  methods.  In  the 
method  of  separation  of  variables  only  geometries  which  coincided  with  the 
separable  coordinate  systems  could  be  treated.  One  often  had  to  drastically 
modify  a practical  problem  in  order  to  fit  it  to  a separable  geometry.  Seme 
of  the  separable  systems  were  extremely  complex  mathematically.  It  is  not 
an  exaggeration  to  sav  that  in  some  cases  years  of  experience  were  re- 


quired for  an  engineer  t.c  effectively  utilize  a particular  coordinate  sys- 
tem. However,  considering  the  techniques  available,  the  only  practical 
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approach  was  to  push  these  limited  methods  to  their  utmost,  and  this  was 
done  most  effectively  during  the  decades  of  the  1930's  and  1940's. 

For  the  practicing  engineer,  the  classical  methods  were  not,  for  the 
most  part,  completely  within  his  grasp;  and  furthermore,  his  problems  were 
much  too  complex  to  be  reduced  to  the  simple  mathematical  forms  which  could 
be  handled  by  the  classical  methods.  As  a result,  he  had  recourse  to  a "bag 
of  tricks"  developed  over  the  years  for  each  particular  area  of  application. 
In  '-■ther  words,  there  was  a severe  "gap"  between  the  simple  applications  of 
the  classical  methods  used  by  researchers  and  the  complex  applications  con- 
fronting the  practicing  engineer. 

During  the  1940's  and  1950' s the  gap  was  closed  somewhat  by  the  advent 
of  the  variational  techniques  which  were  applied  first  to  waveguide  and 
then  to  radiation  p-obiems.  Restrictions  on  the  geometry  were  partially 
removed.  Variational  expressions  could  be  found  for  a wide  variety  of 
quantities  of  interest.  However,  the  evaluation  of  these  expressions  was 
often  quite  difficult  because  of  the  multiple  integrals  involved. 

During  the  decades  of  the  1960's  and  1970' s a new  method  has  emerged 
which  is  capable  of  treating  a wide  variety  of  problems.  This  is  the 
method  of  moments  as  outlined  by  Harrington  [1],  [2].  This  method  has 
proven  especially  useful  because  of  its  wide  applicability.  The  method  is 
conceptually  simple,  although  the  detailed  methodology  may  be  quite  complex. 
In  each  case,  the  defining  integral  equation  is  replaced  with  a matrix  equa- 
tion. The  matrix  equation  may  then  be  solved  by  any  one  of  a number  of 
techniques,  suen  as  matrix  inversion. 

In  the  method  of  moments,  prior  knowledge  of  the  induced  source  distri- 
bution (charge  distribution  for  static  problems  and  current  distribution  for 
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time-varying  problems)  is  not  required,  and,  as  a result,  complex  structures 
can  be  treated  almost  as  easily  as  simple  ones.  The  structure  of  the  antenna 
or  scatterer  may  be  very  complex,  but  the  basic  method  does  not  change. 

Methods  of  matrix  inversion  for  the  solution  of  problems  in  electro- 
magnetic theory  date  back  almost  a century  to  Maxwell  [5]  who  used  such  a 
method  to  calculate  approximate  values  of  capacitance,  and  who  called  the 
technique  "the  method  of  subareas".  In  the  1950's,  Reitan  and  Higgins  [ 6] — [ 10} 
applied  the  same  method  to  a number  of  capacitance  problems.  In  the  early 
1960's,  several  investigators  £11 ] — f 38 ] applied  the  method  to  time-harmonic 
radiation,  scattering  and  waveguide  problems.  Harrington  [1],  [2]  unified 
the  subject  material  and  gave  it  the  name  "the  method  of  moments". 

In  the  late  1960's  and  early  1970' s there  has  been  a tremendous  out- 
pouring of  literature  on  the  subject  of  moment  methods.  The  methods  have 
been  applied  to  an  extremely  wide  variety  of  problems  in  electromagnetic 
theory,  including  thin-wire  radiation  and  scattering  problems,  scattering 
from  dielectrics,  analysis  of  waveguides  of  arbitrary  cross  section,  dis- 
continuity problems  in  waveguides  and  microstrip,  analysis  of  lossy  structures, 
propagation  over  an  Inhomogeneous  earth,  antenna  beam  pattern  synthesis, 
analysis  of  near  fields,  analysis  of  radiation  and  scattering  from  bodies  of 
revolution  and  bodies  of  arbitrary  shape,  etc.  Many  of  these  were  problems 
which  could  not  be  adequately  treated  heretofore. 

The  literature  on  the  method  of  moments  is  already  so  large  as  to  pre- 
clude a comprehensive  bibliography.  References  £39)— [115]  Include  a partial 
bibliography  of  the  work  done  at  Syracuse  University.  Short  courses  in  the 
method  of  moments  have  been  offered  at  the  University  of  Mississippi  [116], [117], 
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the  University  of  Illinois  [118],  Ohio  State  University  [119],  UCLA  [120], 
the  University  of  Naples  [121]  and  the  University  of  Trondheim  [122].  The 
course  notes  from  the  University  of  Illinois  are  now  available  in  book 
form  [123]. 

The  method  of  moments  permits  us  to  model  problems  realistically.  The 
boundary  conditions  are  imposed  directly  and  the  matrix  inversion  procedure 
yields  the  current  distribution  which  fits  the  boundary  conditions  imposed. 

The  principal  advantages  of  the  moment  methods  are  flexibility,  completeness, 
accuracy  and  conceptual  simplicity.  Arbitrary  radiation  shapes  and  arbitrary 
excitation  and  loading  are  permitted.  Currents  are  not  assumed  but  are 
computed  and  all  mutuals  are  taken  into  account.  The  moment  methods  have 
proven  to  be  highly  accurate  for  a wide  variecy  of  problems  associated  with 
radiation  and  scattering  from  thin  wires.  Since  the  methods  are  simple  in 
concept,  they  are  accessible  to  persons  with  a wide  variety  of  backgrounds. 

The  principal  limitation  of  the  method  of  moments  is  associated  with  the 
number  of  unknowns  (n) . Problems  of  moderate  size  (nsf200  or  300)  may  be 
readily  treated,  but  significantly  larger  systems  cannot  be  adequately  treated 
at  present,  due  to  computer  core  and  CPU  requirements. 

Clearly,  the  method  of  moments  hr.3  taken  the  antenna/electromagnetic 
communities  by  storm.  For  the  first  time,  we  can  treat  problems  "as  they 
exist",  with  realistic  modeling  of  actual  structures.  For  the  first  time, 
we  have  a single  method,  capable  of  being  applied  to  the  full  range  of  electro- 
magnetic problems,  and  capable  of  being  simplified  or  extended  according  to 
requirements  of  accuracy. 
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1.2  THE  METHOD  OF  MOMENTS  FOR  PARALLEL  STRAIGHT  THIN-WIRE  RADIATORS 
AND  SCATTERERS 


1.2.1  GENERAL  FORMULATION 

In  this  section  the  method  of  moments  is  introduced  for  the  case  of 
parallel  straight  thin-wire  radiators.  Fig.  1-la  shows  a collection  of  N 
straight,  parallel,  z-directed  antennas,  excited  at  arbitrary  locations  by 
ideal  voltage  generators  at  radian  frequency  u.  First  it  is  desired  to  com- 
pute the  current  distributions  on  the  antennas.  Then,  once  the  unknown 
currents  are  determined,  other  parameters  of  interest,  such  as  impedances, 
far  fields,  near  fields,  etc.  may  be  determined  by  further  calculation.  The 
currents  are,  of  course,  affected  by  interactions  (mutuals)  among  the  various 
radiators  and  these  must  be  taker,  fully  into  account. 

Fig.  1-lb  shows  a dipole  antenna  of  length  £ and  radius  a,  with  surface 


currents  J (longitudinal  current),  J.  (azimuthal  current)  and  J (radial 
z f P 


current).  In  the  general  radiation  or  scattering  problem  the  three  components 


of  current  are  all  present.  Furthermore  J and  J,  varv  with  both  z and  6 and 

z 


Jp  varies  with  p and  <J>.  Foi  thin-wire  radiators  (£  » a ard  a <<  X),  circum- 


ferential (J  ) and  radial  (J  ) currents  may  be  neglected,  the  azimuthal  varia- 
9 P 


tions  of  J may  be  neglected,  and  the  surface  current  J may  be  replaced  with 
a filamentary  model  (Fig.  1-lc).  The  justification  of  these  simplifying 
assumptions  is  discussed  in  Appendix  A.  It  should  be  noted  that  the  boundary 
conditions  are  applied  at  finite  radius  "a"  so  that  the  finite  thickness  of 
the  dipole  may  be  taken  into  account.  The  voltage  source  is  replaced  with  an 
equivalent  current  source  (filament)  in  the  gap  region.  These  foregoing 
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Figure  1-1  Thin-wire  Antennas 


(a)  A collection  of  straight,  parallel,  thin  wire  antennas 

(b)  Fat  dipole  model 

(c)  Filamentary  model 

(d)  Filamentary,  subsectioned  model  of  antennas  in  1-la  .. 
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simplifying  assumptions  are  used  throughout  this  volume. 

The  rigorous  treatment  of  fat  dipoles  is  considerably  more  complex. 

The  isolated  fat  dipole  as  a radiator  or  scatterer  may  be  treated  by  the 
body-of-revolution  methods  outlined  by  Mautz  and  Harrington  [79],  [93],  [96], 
[97],  [IOC].  Recently  chat  method  has  been  extended  to  permit  the  treatment 
of  parallel  fat  dipoles  [24]  as  well  as  the  general  problem  of  parallel  bodies 
of  revolution. 

Thus,  the  thin-wire  antennas  of  Fig.  1-la  may  be  represented  as  filaments, 
as  shown  in  Fig.  1-ld.  The  antennas  are  divided  into  subsections.  The  current 
is  unknown,  but  it  is  assumed  in  this  section  that  it  is  constant  over  each 
electrically  small  subsection.  The  current  on  a typical  subsection  is  repre- 
sented by  i^  where  i^  is  a complex  number.  Fig.  1-2  shows  typica-1  subsections 
A£^  and  AH^,  with  currents  i^  and  1^,  voltages  v^  and  v^ , tangential  electric 
fields  Ezi  and  . midpoints  1(7^, y^z^)  and  j(x^,y^,Zj),  endpoints  i + and 
j +,  and  with  current,  voltage,  and  electric  field  polarities  as  indicated. 

In  general  the  unknown  filamentary  current  may  be  designated  as  i(£) 
where  H is  a coordinate  spanning  the  filaments.  The  current  i(£)  may  be  ex- 
pressed in  terras  of  arbitrary  expansion  functions  as  follows: 


KO  - l ij  fjU) 


(1-1) 


where  i^  are  unknown  complex  numbers  and  f^  are  known  expansion  functions. 
Typical  expansion  functions  and  typical  resultant  current  representations 


are  shown  in  Fig.  1-3.  For  purposes  of  this  section  the  f^  are  chosen  as 
pulse  functions: 


fjOO  “ Pj(£) 


(1-2) 
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where 

p^(£,)  * 1 on  the  jth  subsection 

* 0 elsewhere 

Note  that  the  antennas  need  not  be  center-driven.  Antennas  with  off- 
center  feeds  are  treated  as  readily  as  in  the  centerfed  case.  Also,  sub- 
section lengths  need  not  be  uniform  and  subsection  radii  may  differ.  Note 
that  the  gap  subsections  are  represented  by  current  filaments.  Essentially 
the  subsections  are  small  enough  so  that  concepts  of  network  theory  may  be 
used;  any  source  which  applies  the  proper  voltage  across  the  gap  is  accept- 
able. The  magnitude  of  the  gap  currents  is  t nknown,  but  their  effect  and 
the  effect  of  all  other  currents  if5  to  apply  the  specified  voltages  across 
the  gaps.  We  choose  current  sources  rather  than  voltage  sources  for  simplicity. 

Thus,  we  have  reduced  the  physical  problem  (Fig.  1-la)  to  the  model  shown 
in  Fig.  1-ld  which  consists  of  N filamentary  antennas  divided  into  n total 
subsections  (n  > N) . For  purposes  of  this  section,  the  current  is  assumed 
constant  on  each  subsection.  Later  in  section  1.2.4  the  more  general  repre- 
sentation of  the  current  is  discussed  in  terms  of  arbitrary  expansion  and 
weighting  functions. 

Now,  superposition  and  linearity  mav  be  used  to  express  the  relation- 
ship between  currents  and  voltages.  The  voltage  across  typical  subsection 


A£a  may  be  expressed  as  a sum  of  the  contributions  due  to  all  the  currents; 


v.  » i,z.,  + . . . + i z. 
i 1 il  n in 


£ *jzij  * 0 (on  metallic  subsections) 


J-l 


(1-3) 


* specified  voltages  (on  driven 
subsections,  i.e.  gaps) 
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where  z^  is  the  voltage  across  subsection  i due  to  a unit  current  on  sub- 
section j. 

Equation  (1-3)  13  written  at  eacli  of  the  n subsections,  leading  to  a 
matrix  equation 

[<)  - [z] [i]  (1-4) 

where  [v]  is  an  n-element  column  matrix  with  zero  elements  at  locations  corre- 
sponding to  metallic  subsections  and  non-zero  elements  at  driven  subsections, 
[i]  is  an  n-element  column  matrix  of  the  complex  currents,  [z]  is  an  nxn 
matrix,  called  the  generalized  impedance  matrix.  It  represents  all  the 
interactions  between  the  n subsections  of  the  system  and  may  be  considered 
to  be  the  n-port  open-circuit  impedance  parameters  of  the  system  of  n sub- 
sections. The  impedance  matrix  [z]  depends  only  on  the  geometry  of  the  prob- 
lem (the  locations,  orientations,  and  radii  of  the  various  subsections)  and 
the  frequency  u. 

The  basic  concept  of  the  method  of  moments  is  thus  very  simple  involving 
expansion  of  the  current  in  subsectional  bases,  the  linearity  of  the  system 
of  Fig.  1-ld  and  direct  application  of  the  boundary  conditions.  Many  other 
types  of  problems  in  electromagnetics  also  lead  to  equations  of  the  same 
general  form  (Eq.  (1-4)).  Now  it  remains  to  evaluate  the  typical  element 
z^  of  the  matrix  [z].  This  is  a classical  electromagnetic  problem,  involving 
the  determination  of  fields  due  to  a given  current  distribution.  In  most 
cases,  closed  form  solutions  are  not  available  and  [z]  must  be  determined  by 
approximate  methods. 


The  voltage  across  a subsecrion  mav  be  defined  in  several  different  ways  corre 
sponding  to  different  weighting  functions  [1].  The  electric  field  in  a given 
situation  is,  of  course,  unique.  The  simplest  definition  of  voltage  would  be 
the  product  of  tangential  electric  field  times  subsection  length  (impulsive 
weighting  function). 


11 


1.2.2  DETERMINATION  OF  z FOR  PULSE  EXPANSION  FUNCTIONS  AND  IMPULSIVE 
WEIGHTING  FUNCTIONS. 

In  this  section  methods  for  determining  z^  for  parallel  wire  antennas 
are  outlined  and  two  slightly  different  computational  techniques  are  described. 
Fig.  1-2  shows  the  basic  geometry  consisting  of  two  parallel  subsections  of 
A^  and  tily  Note  that,  from  Eq.  (1-3),  z ^ may  be  represented  as  follows 

vi 

» — (with  all  currents  other  than  i^  equal  to  zero)  (1-5) 

which  corresponds  to  the  usual  definition  of  open-circuit  impedance  in  net- 
work theory.  Thus,  z^  may  be  determined  from  a knowledge  of  the  tangential 
electric  field  E zi  and  the  corresponding  voltage  v^,  at  subsection  Ai.^,  due 
to  a unit  current  on  subsection  A£ ^ . Note  that 


Vi  ” -Ezi  ^*i^ 


(1-6) 


because  of  the  polarities  assumed  in  Fig.  1-2. 

The  problem  is  thus  the  classical  one  of  determining  the  electric  fields 
due  to  a filament  of  current  with  finite  length.  Maxwell's  equations  for  the 
time-hnrmonic  case  are: 


V x £ “ -jwuH 

V x 8 ■ jweE  + J_ 


(1-7) 

(1-8) 


Thr  electric  and  magnetic  fields  may  be  expressed  in  terms  of  the 

magneti-.1  vector  potential  A and  the  scalar  electric  potential  <S  where 

T,  i\  -jkjr-r'j 
Iff  JL(£  > e J 1 

A -W  dv’ 

III.  Atrir  - r’  I 


(1-9) 
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(1-10) 


where  r_  and  jc’  represent  vectors  from  the  origin  to  field  point  (x,y,z)  and 
source  point  (x',y’,z'),  respectively.  The  volume  v' includes  all  sources 
(currents),  k is  the  free  space  wave  number.  is  volume  current  density 
and  p is  volume  charge  density. 

Equations  (1-9)  and  (1-10)  are  obtained  from  Maxwell's  equations  and 
the  following  assumptions: 


B * V x A 
_E  = - jtoA  - V4> 

V4A  = -jweu® 

The  differential  forms  of  eqs.  (1-9)  and  (1-10)  are: 
(V2  + k2)A  - -uJ 
(V2  + k2)$  * -p/e 

The  equation  of  continuity  nay  be  represented  as 
V-J  * -joip 
V-Jg  * -juo 
V*_I  * -juX 


(1-11) 

(1-12) 

(1-13) 

(1-14) 

(1-15) 

(1-16) 

0-17) 

(1-18) 
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where  and  I_  represent  surface  and  line  current  density  and  a and  X represent 
surface  and  line  charge  density.  Note  that  the  subscript  s is  omitted  in  Fig. 
1-lb.  Eq.  (1-17)  applies  to  the  actual  finite-radius  wire  model  and  Eq.  (1-18) 


, ‘ -v‘  «■  * \ <?.«•***#<* 


and  i//(+Az‘/2,  -Az/2)  represents  ^ of  Fig.  1-2  with  subsection  A£.  shifted  as  t 

•*  i 

Az'/2  upwards,  and  subsection  A£^  shifted  as  Az/2  downwards,  etc.  Az',  Az 
represent  increments  used  for  numerical  differentiation.  If  these  increments 
are  chosen  equal  to  corresponding  subsection  lengths  (Az1  * Afc  , Az  ■ Afc^) , 
then  Eq.  (1-29)  becomes 

z.jj  * jtopA^Afcji/*  + [i|»(+»+)  - <K+,“)  - <K“>+)  + tf»(-.-)3  (1-31)  | 

! 

where  ij/(+,-)  represents  t|>  of  Fig.  1-2  with  subsection  A£^  shifted  A£^/2  up- 
wards and  subsection  A£^  shifted  <.£^/2  downwards,  etc.  1 1>  may  be  evaluated  by 
methods  given  in  [1],  [2].  Eq.  (1-31)  agrees  with  Eq.  (22)  of  [3]  derived 
under  different  assumptions:  pulse  expansion  functions  for  both  charges  and 
currents.  Thus  the  two  methods  lead  to  an  identical  formulation  if  the  incre- 
ments Az,  Az'  used  for  evaluating  derivatives  are  subsection  lengths. 

Thus  the  antenna  problem  of  Fig.  la  has  been  reduced  to  a set  of  linear 
equations  which  can  be  solved  by  standard  techniques  for  such  systems.  The 
crux  of  the  problem  is  to  determine  the  generalized  impedance  matrix  [z] 
which  represents  all  the  interactions  among  the  current-carrying  subsections.  j 

i 

Note  that  [z]  depends  only  on  the  geometry  of  the  problem,  [z]  is,  in  effect,  I 

i 

I 

the  open  circuit  impedance  matrix  of  an  n-port  system  consisting  of  the  n j 

t 

subsections.  The  matrix  [z]  is  independent  of  excitation  and  loading.  Any  ! 

! 

subsection  of  the  antenna  system  may  be  opened  up  and  excited  with  a voltage 
generator  or  loaded  with  a complex  impedance,  without  changing  the  matrix  [z]. 

The  effects  of  loading  are  explained  in  more  detail  in  section  two. 
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1.2.3  THE  SCATTERING  PROBLEM 


Consider  the  scattering  problem  shown  in  Fig.  1-4  which  shows  a known 
,inc 


field  E incident  upon  a collection  of  straight,  parallel  thin  wires  which 
have  been  replaced  with  filamentary  models.  The  wires  have  been  divided 
into  n subsections  as  before. 

The  boundary  conditions  are  now  applied  to  the  total  tangential  electric 


field  E . at  the  center  of  each  wire  subsection, 
zi 


Let  E 


£(s)  + £(inc) 


(1-32) 


where  JS^nc^  is  the  known  incident  field  and  JE^  is  the  field  radiated  by 


the  current-carrying  wires. 

Taking  the  z component  of  Eq.  (1-32)  and  specializing  the  field  point 


to  the  center  of  subsection  AA^: 


l 

E , - 
zi 

(s)  (inc) 

zi  “zi 

(1-33) 

l 

R 

Define  incident. 

scattered,  and  total  voltages  as  follows: 

l 

S 

f 

i 

(inc) 

i 

- -E(*nc)  Ail 
zi  i 

(l-34a) 

v^s) 

- -E(^  AA, 
zi  i 

(l-34b) 

Vi 

- -E  . AA, 
zi  i 

(l-34c) 

■ 

and  let  [v^nc^] 

(s) 

, (v  ) , and  [v]  be  the  corresponding  column  matrices  of 

subsection  voltages. 


Xs) 


The  scattered  field  E^  may  be  evaluated  by  Eqs.  (1-29)  - (1-31)  since 
(s) 


the  field  E^  desired  is  that  due  to  radiating  wires,  as  before 

n 

r.(s) 


zi 


" l V 

j-1 


ij 


(1-35) 
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and  writing  Eq.  (1-35)  at  the  center  of  each  subsection  yields  a matrix 
equation 


[v(e)]  - [z][ij 


(1-36) 


The  equivalent  of  Eq.  (1-33)  is 
[v] 


[v(s)l  + tv(1"c)] 


(1-37) 


which  becomes,  for  metallic  scatterers, 
[v]  - 0 - [v(s)]  + [v(inc)] 


(1-38) 


or 


t-v(inC>]  - lz)[l) 


(1-39) 

where  the  typical  element  o.'  [-v^nu^]  is  E^nc^(A£^,).  Eq.  (1-39)  may  be 
solved  by  matrix  inversion  (or  other  techniques)  for  the  unknown  wire  currents  [i] . 

The  scattering  problem  is  thus  reduced  to  the  matrix  equation  (1-39)  which 
is  similar  in  form  to  the  matrix  equation  (1-28)  except  that  the  left  hand  side 
of  (1-28),  the  voltage  excitation  matrix,  has  been  replaced  with  the  negative 
of  the  incident  voltages.  The  reason  for  the  sign  reversal  is,  of  course,  tied 
in  with  the  distinction  between  scattered,  incident,  and  total  fields,  [z] 
includes  all  interactions  among  subsections  of  the  scatterers  and  is  thus  identi- 
cal to  the  [z]  matrix  which  would  be  used  to  analyze  the  radiation  properties 


of  this  system  of  conductors  (the  radiation  properties  of  the  system  of  Fig.  A 


excited  and  loaded  at  one  or  more  subsections) . Note  that  the  incident  or 
impressed  field  may  be  that  of  a plane  wave,  a spherical  wave  or  any  other  type 
of  wave.  The  basic  assumption  is  that  the  source  of  the  wave  is  an  independent 
source  (not  affected  by  the  scattered  radiation  from  the  wires). 
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1.2.4  ARBITRARY  WEIGHTING  AND  EXPANSION  FUNCTIONS 


In  sections  1.2.1  and  1.2.2  the  current  has  been  assumed  constant  ever 


each  subsection  (pulse  expansion  function).  Fig.  1-3  indicates  that  better 


approximations  to  the  current  are  obtained  with  more  sophisticated  expansion 


functions.  Moreover,  point-matching  has  been  used  in  sections  1.2.1  and  1.2.2, 


i.e. , the  boundary  conditions  have  been  applied  at  specific  points  (subsection 


centers).  The  tangential  electric  field  is  zero  at  the  center  of  each  metallic 


subsection,  but  there  is  no  control  over  its  value  at  intermediate  points  be- 


tween subsection  centers.  A more  sophisticated  application  of  the  boundary 


conditions  would  be  to  set  an  average  or  weighted  average  of  the  tangential 
field  equal  to  zero  at  metallic  subsections  in  radiation  problems. 

Consider  the  electric  field  E due  to  currents  on  the  wires  of  Fig.  1-1 
as  given  by  Eq.  (1-24): 


±_  (ii_ + p2  f y*'*. e 

ue  \ 2 + k 5 I 4ir  |r , - r’ 


zi  jwe 


wires 


- 


let  I(z')  ■ £ I^f^(z),  as  in  Eq.  (1-1)  where  the  f^  are  arbitrary  expansion 
functions,  and  let  w(z)  be  an  arbitrary  (weighting)  function.  Substituting  in 


Eq.  (1-24): 


2 

E m J X - (- — + k^)  I 


-)k}r.  -r'  I 

2,  f lli!2  e 

) 1 *»ltj  - r'| 

k a 


(1-40) 


Multiplying  both  sides  of  Eq.  (1-40)  by  w(z)  and  integrating  over  subsection 


^Eziw(z)dz-[  I3  J jSI(^2+k2)dz  | -1~-[7  - r»l  d2>  (1~41) 


-JklM'l 


j-1  A£ 


' l';S:f^Si'S/-'-i‘'':  r Jrwe^A^Tidij'? 


Interchanging  the  order  of  differentiation  and  integration: 

-Jk|r.-r’ 


f r f f 1 a2  2 fj(z,>w(z>  e 

J E.iwWdl  - £ TJ  J J jST  (-J  + k > — 4»|V,  -"£•' 

At.  1-1  At.  At.  i 


.1-1  ^Atj  At  j 


dz'dz 




(1-42) 


‘ij 


now  Eziw(z)dz  is  defined  as  voltage  v^,  the  double  integral  on  the  right 

hand  side  of  Eq.  (1-42)  is  defined  as  the  generalized  impedance  z^,  and 

Eq.  (1-42)  becomes 
n 

vi*  1 Vj 
J-l 

which  in  turn  becomes 

M - [z][i]  (1-43) 

as  i runs  through  the  indices  1 - n. 

Note  that  the  voltage  definition  has  been  generalized  and  now  represents 
a weighted  Integral  of  the  electric  field.  If  w(z)  is  an  impulse  and  f^(z') 
is  a pulse  function,  then  Eq.  (1-42)  reduces  to  Eq.  (1-26).  A furthar  generali- 
zation of  Eq.  (1-42)  may  be  made  by  replacing  the  Integrals  with  inner  products  [1] . 

1.3  THE  METHOD  OF  MOMENTS  FC  i NON-PARALLEL  THIN  WIRES 

The  development  of  sections  1.2.1  and  1.2.2  may  be  extended  to  non-parallel 

2 o 2 

(skewed)  wires  (Fig.  1-5).  This  involves  replacement  of  (3  /3z^  + k ) with  the 
operator  (7  x V x ) and  the  use  of  various  trigonometric  relationships.  The 
development  is  laborious  but  straightforward.  Details  of  the  development  for 
pulse  expansion  functions  are  given  in  the  literature  [110]. 

The  result  for  z^j  is 

*ij  " “ Jmc  cos  0ij*  " W+»+)  **  *(-.♦>  ~ *<+»->  + *(“.->11  (l-44)j 
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II  - LOADING  AND  EXCITATION  IN  THE  METHOD  OF  MOMENTS 


In  this  ciiaptcr  the  effects  of  excitation  and  loading  of  subsections 
are  considered.  Several  different  methods  for  treating  loads  are  outlined. 

2.1  SUBSECTION  (n-PORT)  PARAMETERS 

Consider  a collection  of  thin-wire  antennas  such  as  that  shown  in  Fig.  1-1. 
The  antennas  are  subsectioned  and  it  is  assumed  that  the  subsections  are  small 

A £ • 

electrically,  i.e.  « 1,  where  A l is  subsection  length.  It  is  assumed  that 
there  are  n total  subsections,  as  in  Fig.  1-ld,  and  that  N subsections  are 
excited  and/or  loaded. 

Figure  2-1  shows  a typical  subsection  Ai4  driven  with  an  applied  volt- 

X 

age  and  loaded  with  a lumped  impedance  element  . v’  and  z^  represent 


the  Thevenin's  equivalent  circuit  of  the  devices  and  loads  attached  to  sub- 
section Al^.  Since  the  subsections  are  small  electrically,  the  response  of 
the  system  of  antennas  to  the  excitation  and  loading  does  not  depend  on  the 
fine  details  of  the  generators  and  load  construction,  but  only  on  the  values 
of  the  applied  voltage  v^  and  load  z^.  As  far  as  terminal  behavior  is  con- 
cerned, the  electromagnetic  problem  becomes  a network  problem.  The  applica- 


tion of  boundary  conditions  leads  to  the  matrix  equation 
M - [z][i] 


(2-1) 


and  the  equations  of  constraint 
[vj  - (v'j  - [z1mj 


(2-2) 
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Combining  (1)  and  (2): 
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or  tv']  - [[z]  + [z^J]  [i]  . [ZjJHi] 
where  [z^J  - [z]  + [z^] 

Solving  for  the  currents: 

[i]  - I^J'V]  - [y’Hv'] 


where 
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and  then  the  total  voltages  (v]  nay  be  determined  by  Eq.  (2-2).  Note  that 
the  effect  of  loading  is  merely  to  replace  matrix  [z]  with  [[z]  4-  [z^]]  and 
[v]  with  [v*]. 

In  Eq.  (2-2),  zeros  occur  on  the  left  hand  side  if  the  total  voltage 
across  a subsection  is  zero.  In  Eq.  (2-3),  zeros  occur  on  the  left  hand 
side  only  if  the  applied  voltage  is  zero.  Thus,  if  a subsection  is  loaded 
but  not  driven,  a zero  would  occur  in  the  appropriate  element  of  [v'].  All 
problems  involving  loading  and  excitation  may  be  treated  by  the  above  formu- 
lation. However,  in  some  cases  a simpler  matrix  equation  may  be  used.  These 
special  cases  are  described  in  sections  2.2  and  2.3. 


2.2  N-PORT  PARAMETERS  OF  THE  DRIVEN  OR  LOADED  SUBSECTIONS 

The  addition  or  changing  of  a load  changes  all  the  currents  and,  as  a 
result,  the  near  snd  far  fields.  In  order  to  determine  these  quantities, 
then,  we  need  to  invert  a new  nxn  matrix  as  in  Eq.  (2-4).  If,  however,  we 
are  concerned  only  with  behavior  at  the  N driven  and/or  loaded  ports,  such 
as  port  voltages,  port  currents,  coupling  from  one  port  to  another,  etc.. 


then  the  inversion  of  an  N by  N matrix  is  sufficient. 
Let  [y]  » [z]  1 


(2-5) 


From  the  larger  n x n [y]  matrix,  the  terms  corresponding  to  the  N driven 
ports  are  selected  to  form  an  N-port  short  circuit  admittance  matrix  denoted 
by  [Y]  (see  (3),  [54],  [129]).  The  N-Port  open  circuit  impedance  matrix, 
denoted  by  [Z],  is  the  inverae  of  [Y]. 


[Z]  - [Yf A 

and  the  voltages  and  currents  at  the  N driven  and/or  loaded  ports  are 


(2-6) 


related  by 


[V]  - [Z][I] 


(2-7) 


The  total  voltages  may  than  be  obtained  from  Eq.  (2-7).  Comparing  Eqs.  (2-4) 
and  (2-10),  it  is  noted  that 
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for  chose  ports  i,j  which  are  loaded  and/or  driven.  Thus  the  matrix  [Y^] 
may  be  formed  by  selecting  appropriate  elements  of  [y ],  as  is  the  case 

X* 

with  [Y]  and  [yj. 


2.3  THE  DUAL  ANALYSIS  FOR  LOADING 

Fig.  2-2  shows  the  Norton's  equivalent  circuit  of  the  devices  and  loads 
attached  to  subsection  Ai^  where  1^  represents  short  circuit  current  and  Y£^ 
represents  the  admittance  seen  when  the  current  generator  is  deactivated. 

Now  we  may  write  the  following  matrix  relationships  between  currents 
and  voltages  at  the  n subsections: 

[i]  - [i*i  - ly£][v]  (2- 


ti]  - [yJM 


where 


vY*u 


v'Jhn 


Then,  transposing  terms  in  (2-13) 


and 


U’]  - tlvl  + [y^l  [v]  = [v'^l  [v] 


(2-16) 


[v]  - tvjjji-1!!']  - I^itri 


(2-17) 
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where  the  subscript  D indicates  a dual  quantity.  Once  the  specified  driving 
voltages  (v’J  are  known,  then  [ i * ] mav  he  determined  from  (2-14),  [v]  may 
be  determined  from  (2-17),  and  [i]  mav  then  be  determined  from  (2-13).  Note 
that  [y^DJ,  [z^J  are  not  necessarily  enual  to  [v^],  [z^]. 

Similarly,  the  following  dual  relationships  mav  readily  be  developed 
for  the  N-Port  analysis: 


and 

where 


[I’l  - ICY]  + [Y£]](V] * [YJd]  [V] 

m - ir^fhv)  - tz;Dni’] 

[I]  - [I*]  - [y£](v] 
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In  chapter  two,  the  effects  of  loading  for  a system  of  n subsections 

have  been  considered.  Mathematically,  the  effect  is  merely  to  replace  the 
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generalized  (n  x n)  impedance  matrix  [z]  with  fz]  + [z^]  where  fz^]  is 
a diagonal  load  matrix,  and  to  replace  [v],  the  total  '-oitaga  excitation 
matrix,  with  [v*],  the  applied  voltage  excitation  matrix.  The  terminal 
behavior  (terminal  voltages  and  currents)  at  driven  and/or  loaded  ports 
may  be  described  by  an  N x N matrix  equation  where  N is  the  total  number 
of  such  ports.  Finally,  the  dual  formulation  for  loading  is  described. 
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LII  - SUMMARY  01’  METHODS  FOR  THE  DKTERMT NATION  OF  AUXILIARY  OUANTIIIES 


3.1  INTRODUCTION 

In  sections  one  and  two,  the  matrix  equation  (1-28)  and  its  related 
forms  (1-39),  (2-3),  (2-16)  have  been  considered  as  the  basic  formulation 
of  the  radiation  or  scattering  problem,  from  which  the  currents  may  be 
determined.  In  this  short  section,  the  determination  of  other  quantities 
of  interest  is  described  briefly  and  references  are  given  for  more 
complete  descriptions. 


3.2  IMPEDANCE  AND  ADMITTANCE  QUANTITIES 

Once  the  currents  and  voltages  are  known  for  a given  system,  input 
impedances  at  a driven  port  are  given  simplv  bv  the  ratio  of  voltage  to 
current  at  that  port.  The  input  impedance  at  anv  subsection  of  the  system 
may  thus  be  determined  by  opening  up  the  subsection  and  exciting  it  with 
a voltage  generator.  The  input  impedance  mav  also  be  determined  directly 
from  [z]  and  [z^]. 

Consider  an  antenna  system  such  as  that  of  Fig.  Id  consisting  of  n 
total  subsections.  Suppose  that  N subsections  are  driven  and/or  loaded. 
The  open-circuit  impedance  matrix  and  the  short  circuit  admittance  matrix 
of  the  N-port  svstem  can  be  readily  obtained  from  [ z ] . First  [y]  is  ob- 
tained by  inverting  [ z ] . Then  the  elements  are  selected,  where  both 
i and  j are  members  of  the  N-port  svstem  (driven  and/or  loaded  ports),  to 
form  an  N x N matrix  [Y]  (See  [3],  [!>4],  [129]).  This  matrix  [Y]  is  the 
short-circuit  admittance  matrix  of  the  N-port  system,  and  its  inverse  [Z] 


33 


&*** 


g*~t *±-f- & •- 


is  the  open-circuit  impedance  matrix.  If  N is  much  smaller  than  n,  it 
may  be  convenient  to  utilize  [Y]  and  [Z]  to  establish  the  terminal  be- 
havior for  different  excitations.  Anv  of  the  other  common  network 
parameters,  such  as  the  scattering  matrix,  AHCI)  parameters,  etc.,  may 
also,  of  course,  be  obtained  from  [Y],  [Z]. 


3.3  NEAR  AND  FAR  FIELDS 

Once  the  currents  are  known,  the  near  and  far  fields  may  be  obtained 
bv  further  matrix  operations.  The  near  fields  are,  of  course,  considerably 
more  complicated  than  the  far  fields.  At  a point  in  the  near  field,  there 
are  ir-  general  three  components  cf  E and  three  components  of  H.  The 
electric  and  magnetic  vectors  at  a point  trace  out  independent  polariza- 
tion ellipses  in  time  [130].  However,  the  basic  information  is  available 
from  the  matrix  [zj.  Since  z^  is  proportional  to  the  electric  field  (in 
a given  direction)  at  field  point  i,  due  to  the  current  on  subsection  j, 
one  need  only  sum  over  all  subsections  to  obtain  one  component  of  .E  at  near 
field  point  i.  The  other  components  of  E are  obtained  in  a like  manner 
and  the  three  components  of  ji  are  obtained  bv  numerical  differentiation  of 
the  $ functions.  Details  of  the  near  field  procedure  are  given  in  [71], 
[73],  [83],  [85],  [89],  [981,  [104],  [105],  [110],  [125]-[128],  [13Q]-[134]. 

The  far  field  computation  represents  merely  a specialization  of 
z^  to  the  far  field  case.  The  far  fields  are  obtained  simplv  by  summing 
the  contributions  of  all  subsections  (See  [3]  for  details).  Directive 
gain  may  be  obtained  by  an  integration  of  the  far  field  beam  pattern. 
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3.3.1  Radiation  Hazards 

Radiation  hazards  may  be  evaluated  at  a point  In  the  near  field  by 
considering  all  complex  components  of  the  near  fields  [130].  The  computer 
program  described  In  [73]  has  a sumoutine  which  computes  radiation  hazards 
in  the  near  field  of  a collection  of  thin-wire  antennas  or  scatterers. 


3.3.2  Near  Field  Ellipses 

Once  the  six  complex  components  of  _E  and  H are  known  at  a near  field 
point,  the  description  of  the  near  field  ellipses  may  readilv  be  obtained 
[130].  Complete  information  on  the  near  field  E and  11  ellipses  (ellipti- 
city,  plane  of  the  ellipse,  orientation  of  the  major  axis)  is  provided  by 
a subroutine  of  the  computer  program  described  ir.  [73]. 


3.4  BEAM  PATTERN  SYNTHESIS 

Beam  pattern  synthesis  can  be  utilized  in  conjunction  with  matrix 
methods.  Given  an  array  of  N antennas  with  N feed  points,  one  may  specify 
N points  of  the  far  field  beam  pattern  and  then  determine  the  required  N 
voltage  excitations  at  the  feed  points.  One  may  also  specify  the  beam 
pattern  at  more  than  N points  and  find  the  N voltage  excitations  for  least 
square  error  with  respect  ot  the  specified  pattern.  The  beam  pattern  may 
also  be  synthesized  subject  to  other  constraints.  Details  on  synthesis 
procedures  are  given  in  [46],  [87],  [151]. 

3.5  OPTIMIZATION  OF  RADIATION  CHARACTERISTICS 

Many  of  the  typical  performance  characteristics  of  an  antenna  system, 

such  as  gain,  efficiency  and  qualitv  factor,  mav  be  optimized  with  or 

without  constraints.  For  instance,  given  a system  of  N antennas,  with  N 
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feed  points,  one  may  readily  obtain  the  value  of  the  maximum  gain  in  a 
given  direction,  and  the  required  voltage  excitations  (at  the  N feed  points) 
which  produce  chat  gain.  Details  of  the  optimization  procedures  are  given 
in  [1],  [3],  [54],  [108]. 


3.6  APERTURE  COUPLING 

Electromagnetic  coupling  through  apertures  in  plane  screens  may  be 
evaluated  by  the  method  of  moments.  Babinet’s  principle  is  employed  tc 
transform  the  problem  to  one  of  scattering  by  the  complementary  disk.  For 
infinitely  long  slots,  the  problem  simplifies  and  mav  be  treated  accurately 
for  slots  of  considerable  width  (<  10X).  The  surface  of  the  complementary 
disk  is  modeled  exactly  in  this  case.  For  apertures  of  arbitrary  shape, 
and  of  moderate  area  ( <*  one  square  wavelength) , the  complementary  disk 
may  be  modeled  approximately  as  a wire  grid.  Then  the  usual  thin-wire 
methods  may  be  applied  to  compute  such  quantities  as  aperture  fields, 
fields  radiated  beyond  the  aperture,  transmission  coefficient,  etc. 

3.7  ANTENNA-TO-ANTENNA  COUPLING 

Antenna-to-antenna  coupling  mav,  of  course,  be  evaluated  directly 
once  all  loads  have  been  taken  into  account  and  all  currents  and  voltages 
on  the  wire  structures  are  known.  The  power  out  of  a driven  section  and 
the  power  absorbed  by  a loaded  subsection  may  be  evaluated  directly  in 
terms  of  the  subsection  current  and  voltage. 

3.8  ELECTROSTATIC  APPLICATIONS 

The  method  of  moments  may  be  applied  to  electrostatic  problems.  The 
details  of  the  development  for  electrostatics  are  given  in  [39],  [40].  The 
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excitations  are  the  static  voltages  applied  to  a system  and  the  response 
is,  for  instance,  the  unknown  induced  charge  distribution.  The  charge- 
bearing surfaces  are  divided  into  subsections  and  appropriate  expansion 
functions  for  the  unknown  charge  distributions  are  assumed.  A matrix 
equation  of  the  form 


£VJ  = [D]  [o ] 


(3-1) 


results,  where  [V]  is  a column  matrix  representing  the  voltage  excitation, 
[oj  is  a column  matrix  representing  the  unknown  induced  charges,  and  [D] 
is  a square  matrix  representing  the  interactions  among  subsections.  D^, 
for  instance,  may  represent  the  potential  at  subsection  i due  to  a unit 
charge  on  subsection  j.  Various  electrostatic  calculations  may  be  carried 
out,  such  as  the  determination  of  capacitance  and  resistance,  and  the 
evaluation  of  electric  fields  and  potentials  [49],  [39],  [40].  The 
method  of  moments  has  been  used  extensively  in  connection  with  the  evalua- 
tion of  the  static  properties  of  microstrip  transmission  lines. 
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IV  - SOME  APPLICATIONS  OF  THE  METHOD  OF  MOMENTS  FOR  THIN  MIRES 

The  method  of  moments  may  be  applied  to  a wide  variety  of  antenna 
and  scattering  problems.  Examples  of  computed  and  experimental  data  for 
representative  problems  are  presented  in  this  section.  Most  of  the  ex- 
amples presented  represent  phenomena  which  could  not  in  the  past  be 
adequately  treated. 


4.1  DIPOLE- TO-DIPOLE  COUPLING 

One  may  compute  tha  coupling  between  any  two  thin-wire  antennas,  using 

the  method  of  moments  and  the  formulation  of  section  two.  Coupling  is 

defined  here  as  the  ratio  of  output  power  to  input  power. 

K.  Siarkiewicz  has  made  an  extensive  theoretical  and  experimental 

study  of  dipcle-to-dipole  coupling  [135]-[137].  Figs.  4-1  to  4-3  show 

* 

some  typical  results  of  this  study.  Fig.  4-1  shows  computed  and  experimental 

data  for  the  coupling  between  two  resonant  parallel  dipoles  oriented  end- 

to-end  as  shown  in  Fig.  4-1.  One  dipole  is  centerfed  and  the  other  is 

center loaded  with  50H.  Dipole  radius  is  0.0119A  and  dipole  length  is  0.48A. 

Note  that  the  loaded  antenna  is  located  at  a beam-pattern  null  of  the  driven 

antenna.  Only  the  radial  component  E^  of  the  electric  field,  which  falls 
2 

off  as  1/r  , contributes  to  coupling.  Accordingly,  the  coupling  approaches 
a 12  db/octave  asymptote.  The  experimental  data  shows  reasonably  good 
agreement  with  the  predicted  results.  There  is  a small  periodic  experi- 
mental error  present,  which  arises  from  ground  reflections,  as  explained 
later.  The  magnitude  of  the  error  increases  with  spacing  d. 

Fig.  4-2  shows  computed  data  for  dipole-to-dipole  coupling  for  A/4 
dipoles  arranged  (a)  end-to-end  and  (b)  broadside.  Dipole  radius  is 


The  ratio  of  output  to  incident  power  is  plotted  in  Figs.  4-1  to  4-3.  This 
differs  from  coupling  if  the  input  reflection  coefficient  is  nonzero. 


Figure  4-1  Coupling  between  resonant  dipoles  oriented  end-to-end. 


^ v -v-  '■*'’¥&$**"  * '"  "'******'*?' 

£'  k ^ 


:«ASPB08Cl^i»ww«  y*^i&&*ttBHts**Z»f*fB^^ 


0.01055X.  The  end-to-end  case  approaches  a 12  db/octave  asymptote,  whereas 
the  broadside  case  approaches  a 6 db/octave  asymptote.  Note  that  for  very 
small  values  of  D/A,  the  coupling  in  the  end-to-end  case  actually  exceeds 
that  of  the  broadside  case. 

Coupling  in  the  presence  of  parasitic  wires  may  also  be  accurately 
predicted  by  moment  methods.  Fig.  4-3  shows  two  resonant-length  dipoles, 
arranged  end-to-end,  with  a resonant  parasitic  element  located  nearby.  In 
effect,  there  are  two  coupling  paths,  one  direct  and  one  by  way  of  the 
parasitic  element.  Both  the  currents  on  the  driven  element  and  those  on 
the  parasite  contribute  to  coupling.  As  spacing  D is  varied,  the  two  contri- 
butions, which  are  nearly  equal  in  magnitude,  add  for  some  spacings  and 
partially  cancel  for  others,  giving  rise  to  the  oscillatory  coupling  effect 
shown  in  Fig.  4-3.  The  measured  data  shows  excellent  agreement,  except 
near  the  nulls.  Note  that  both  the  positions  and  values  of  the  peaks  in 
coupling  are  accurately  predicted,  despite  the  rapid  oscillation  as  a 
function  of  D/A.  A small  experimental  error  due  to  ground  reflections  is 
believed  to  be  responsible  for  the  departure  near  the  nulls.  One  would 
not,  of  course,  expect  good  agreement  near  the  nulls  if  there  were  any 
reflections  present.  The  experimental  error  caused  by  ground  reflections 
increases  with  increasing  D/A. 

The  experiments  for  Figs.  4-1  and  4-3  were  performed  by  elevating  the 

dipoles  about  15A  above  a concrete  floor.  Small  reflections  from  the  floor 

produce  interference  at  tlje  deep  nulls  of  Fig.  4-3  and  also  give  rise  to 

the  small  periodic  experimental  error  shown  in  Fig.  4-1.  As  separation 

of  the  dipoles  increases,  ,the  level  of  the  reflected  signal  becomes  more 
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significant  with  respect  to  the  direct  signal.  Thus  experimental  errors 
increase  with  dipole  separation,  as  indicat'  in  Figs.  4-1  and  4-3. 

Siarkiewicz  has  shown  [ 135]  that,  for  parallel  dipoles  on  a common 
ground  plane,  coupling  may  be  accurately  predicted  by  a simple  empirical 
formula.  He  has  also  shown  that,  for  electrically  small  dipoles,  coupling 
cannot  be  accurately  measured  (special  corrections  which  utilize  theoreti- 
cal data  are  required) . More  complete  information  on  coupling  computations 
is  included  in  references  (135]- [137]  and  in  a forthcoming  volume  of  the 
Method  of  Moments  Applications  Series. 


4.2  NEAR-FIELD  PREDICTION 

The  study  of  near  fields  has  bean  of  interest  to  antenna  engineers  for 
many  years.  It  has  been  a difficult  study  because  of  the  generally  complex 
nature  of  the  near-field  structure.  Even  for  thin-wire  antennas,  there  are 
only  a few  specific  cases  that  have  been  studied  in  detail.  However,  in- 
terest in  such  problems  is  high  because  of  the  number  of  relatively  high- 
power  wire  and  rod  antennas  and  high-gain  arrays  in  use  and  because  of  con- 
cern over  possible  associated  radiation  hazards  and  compatibility  problems. 

The  near  fields  of  individual  dipoles  and  linear  dipole  arrays  have 

been  studied  in  some  detail  and  some  of  the  results  are  '-'nown  in  F-*gs.  4-4 

to  4-9.  Fig.  4-4  shows  the  near  electric  and  magnetic  fields  of  a center- 

fed  halfwave  dipole  of  a radius  a “ 0.0CrX.  Voltage  excitation  is  one  volt 

and  frequency  is  100  MHz.  Fig.  4-4a  shows  the  radial  component  E of  the 

P 

electric  field,  as  computed  by  two  different  techniques  of  the  method  of 
moments.  Method  1 uses  pulse  expansion  functions  and  method  2 uses  tri- 
angular expansion  functions.  The  computer  programs  are  described  respectively 
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Figure  4-4a  Near  fields  Ep  of  a halfwave  dipcle  (a/A  = 0.005,  A£/A 
V « I volt,  A = 3m). 
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in  references  [73]*  (69].  Also  shown  in  Fig.  4-4a  is  an  approximate  calcu- 
lation of  near  fields  from  current  distributions,  based  on  a generalized 
Gauss'  law  formulation  [71].  Note  the  excellent  agreement  among  the  three 
different  methods,  even  at  distances  as  close  as  0.01A  (for  o/a  “2).  E 


has  peaks  near  the  gaps  and  near  the  ends  of  the  dipole.  In  the  central 


region  (z;5s  0.125A) , has  a (1/p)  dependence  on  p.  Methods  1 and  2 


diverge  slightly  near  the  ends  of  uhe  dipole  and  this  is  expected  because 
of  the  difference  in  expansion  functions.  Fig.  4-4b  shows  the  tangential 
component  E^  of  the  electric  field  as  computed  by  methods  1,  2.  E^  has  ap- 
proximately a £n(p/a)  dependence  in  the  central  region  (z  * 0.125A).  Fig. 
4-4c  shows  the  magnetic  field  component  H^,  as  compared  with  an  approximate 
calculation  of  near  fields  from  current  distribution,  based  on  a generalized 
Ampere's  law  formulation  [71].  has  approximately  a 1/p  dependence.  The 
data  shown  In  Fig.  4-4  have  also  been  compared  with  other  results  in  the 
literature.  The  computed  data  on  near  fields  E , E , H-  show  very  close 

P 2 P 


agreement  with  Harrison  et  al, [138]  for  p - 0.05A.  The  data  for  E lies 

z 


between  those  of  King  and  Wu  [139]- [141]  and  Galejs  [142]  for  p < 0.05A. 

Fig.  4-5  shows  a comparison  of  dipole  neat-field  computation  with  ex- 
perimental results.  A third  technique  of  the  method  of  moments,  that  of 
ref.  [127],  with  sinusoidal  expansion  functions,  is  added.  Figs.  4-5a,  b 


show  data  for  E , E respectively  at  p ■ 0.03A.  Even  at  this  close 
P 2 


spacing,  the  agreement  among  computational  methods  and  between  theory  and 
experiment  is  excellent. 

Fig.  4-6  shows  the  near  electric  field  of  an  eight-element  linear  array 


of  centerfed  halfwave  dipoles  with  halfwave  spacing.  Fig.  4-6a  shows  the 
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Figure  4-5a  Near  fields  E of  a resonant  dipole  - theory  and  experiment 

P 

(a/A  - 0.00503,  £/A  - 0.04484,  A£/A  - 0.02242,  A - 3m). 
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Figure  4-6b  Near  electric  field 
d/A  = 0.5,  z=0)~ 

of  an  eight-element  array  of  dipoles  (£/X  *»  0.5, 
Polar  coordinate  plot,  endfire  excitation. 
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near  electric  field,  in  the  principal  H-plane,  for  uniform  excitation.  Note 
that  the  patterns  smoothly  and  gradually  approach  the  far  field  pattern  with 
the  peaks  and  nulls  forming  gradually  as  r (radial  distance  to  the  field 
point)  increases.  The  pattern  for  r ■ 40A  is  extremely  clos':  to  that  of 
the  far  field.  Since  mutuals  are  taken  into  account,  there  is  some  null 
fill-in  and  the  far  field  pattern  differs  slightly  from  the  classical 
pattern  of  an  array  of  Hertzian  dipoles.  Fig.  4-6b  shows  the  near  electric 
field  of  the  eight  element  array,  in  the  principal  H-plane,  for  endfire 
excitation.  Note  that  the  main  beam  broadens  as  radial  distance  r increases, 
in  contrast  to  the  broadside  case  (Fig.  4-6a)  which  shows  main  beam  narrowing 
as  r Increases. 

For  points  in  the  "very  near"  field,  the  electric  and  magnetic  fie.'ds 
are  strongly  effected  by  the  mutual  interactions  among  dipoles  [71].  Fig. 

4-7a  shows  a plot  of  electric  field  strength  across  the  face  of  the  eight- 
element  array,  at  distances  between  .05  and  .251.  Fig.  4-7b  shows  a corre- 
sponding plot  of  the  electric  field  strength  of  an  eight-element  array  of 
Hertzian  dipoles.  Apparently  the  mutuals  have  the  effect  of  "smoothing  out" 
the  variations  of  electric  field  in  the  "very  near"  region.  The  r or responding 
plot  for  the  Hertzian  dipoles  in  the  region  r - 3 - 40A  differs  only  slightly 
from  that  of  Fig.  4-6a. 

The  radiation  hazards  in  the  near  fields  may  be  computed  [71],  [130]  once 
all  the  components  of  the  electric  field  are  known.  Radiation  hazard  contours 
may  readily  be  drawn  from  the  printed  output  of  a so-called  "snapshot"  routine 
[73]  which  prints  out  radiation  hazard  values  over  a given  rectangle  in 
space.  Fig.  4-8  shows  an  example  of  such  a plot  in  the  principal  H-plane  (z 
■ 0)  of  the  uniformly  excited  eight-element  array.  Each  dipole  is  excited 
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Figure  4-7a  Near  electric  field  of  a uniformly  excited  array  of  dipoles  (d/X  * 0.5, 
z = 0)  - Rectangular  plot,  halfwave  dipoles. 
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Figure  4-8  Radiation  hazard  contours  in  the  principal  H-plane  of  a uniformly 

excited  eight-element  linear  array  of  dipoles  (v  = 1 volt,  l/X  «•  0.5, 
d/X  » 0.5,  z » 0) . 
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with  one  volt.  The  assumed  radiation  hazard  criterion  is  10  mw/cm  or 


194  volts/m.  Fig.  4-8  shows  decibel  levels  relative  to  that  criterion. 


Note  that  the  partially  formed  peaks  and  nulls  of  the  near  field  beam  patterns 


(Fig.  4-6a)  show  up  clearly  in  the  plot  of  Fig.  4-8,  forming  progressively 


(as  r increases)  at  angles  corresponding  to  peaks  and  nulls  of  the  far  field 


beam  pattern. 


Fig.  4-9  shows  a comparison  of  theory  and  experiment  for  the  near-field 


component  E in  the  principal  H-plane  of  a two-dipole  endfire  array  of  resonant 
z 


dipoles  with  halfwave  spacing.  The  agreement  between  theory  and  experiment  is 


good  except  at  positions  very  close  to  the  ends  of  the  dipole. 


More  complete  information  on  near-field  problems  is  given  in  references 


£71] , [73],  [ 130]  — [ 133 ] , and  in  volumes  If 2,  If 3,  If 5 of  the  RADC  Method  of 


Moments  Applications  Series. 


4.3  TRANSMISSION  '.IKE  INTERFERENCE  WITH  DIPOLE  PERFORMANCE 


Fig.  4-10  shows  a transmission  line  of  length  wire  radius  "a1 


wire  center  separation  d,  connected  to  a dipole  of  length  £ and  wire  radius 


"a".  The  transmission  line  interferes  with  the  dipole;  it  contributes  a cross 


polarized  component  E^  as  well  as  interfereing  with  the  normal  dipole  pattern. 


The  characteristics  of  such  thin-wire,  finite-gap  A/2  and  A/4  dipoles 


fed  with  a length  of  two-wire  transmission  line  have  been  studied  using 


matrix  methods  [75],  [113].  Transmission  line  length  and  gap  size  have  been 


varied  obtaining  the  complete  solution  (current  distributions,  impedance  and 


far  field  beam  patterns)  in  each  case.  Significant  departures  from  the  ideal 


dipole  characteristics  are  found  even  for  electrically  small  gap  sizes  and 


transmission  line  lengths.  When  the  departrres  are  small,  they  can  be  described 
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Figure  4-9a  Near  electric  field  of  a two-element  endfire  array  of  resonant 
dipoles  U/A  = 0.4484,  d/A  = 0.5,  a/A  » 0.00503,  A£/A  = 0.02242, 
A ■ 3m,  V = 1 volt)^  y/A  * 0.05. 
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Figure  4-10  Dipole  and  Transmission  line. 
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adequately  by  a cross  polarization  ratio  and  by  null  fill-in.  These  effects 


are  minimized  by  choosing  transmission  line  length  £t  equal  to  nX,  where  n 


is  an  integer.  Departures  from  ideal  dipole  performance  are  significantly 


increased  as  dipole  length  i decreases,  gap  length  d increases,  transmission 


line  length  fc  increases  and  current  imbalance  increases.  The  cross  polari- 


zation ratio  (defined  as  |Eft  | / j |)  is  minimized  by  choosing  line 

v max  Tj  max 


length  ij.  equal  to  nA  where  n is  an  integer.  Fig.  4-11  shows  computed  data 


for  the  cross  polarization  ratio  and  Fig.  4-12  shows  a comparison  with 


experimental  data.  Note  the  deep  null  in  the  cross  polarization  ratio  for 


=*  1.0.  More  complete  information  on  transmission  line  interference 


with  dipole  performance  is  given  in  references  {75],  [144]  and  in  volume 


// 4 of  the  RADC  Method  of  Moments  Applications  Series  [113]. 


4.4  MICROSTRIP  APPLICATIONS 


Static  microstrip  problems  are  ideally  suited  to  treatment  by  the 


method  of  moments  and  the  results  mav  be  used  for  "quasi-TEM"  analysis  of 


microstrip.  The  static  capacitance  matrix  can  be  computed  by  three  methods 


associated  with  the  method  of  moments.  In  the  first,  unknown  free  and 


bound  charges  are  assumed  on  all  interfaces  and  simple  free-space  Green's 


functions  are  employed.  The  basic  matrix  equation  includes  boundary  condi- 


tions on  potential  over  each  conductor  and  boundary  conditions  on  continuity 


of  the  normal  displacement  vector  ovar  each  dielectric  interface.  Complex 


problems  involving  multiple  dielectrics,  conductors,  and  ground  planes  can 


be  treated  bv  this  method.  However,  the  number  of  unknowns  in  the  matrix 


equations  may  be  large.  This  method  has  been  described  by  Adams  and  Mautz 


[49],  [40]  and  was  implemented  in  a computer  program  by  Priebe  and  Kyle  [145] 


CROSS  POLARIZATION  RATIO  (db) 


CROSS  POLARIZATION  RATIO  (db) 


x v *»*«»•»! 


for  arbitrary  raicrostrip  geometries.  Fig.  4-13  shows  the  capacitance  of 
a dielectric-loaded  parallel-plate  computed  by  this  method.  Results  for 
■ 1 are  compared  with  those  of  Harrington  [1]  and  Reitan  [10]. 

In  the  second  method,  a Green's  function  for  the  basic  microstrip 
geometry  is  obtained  by  images,  leading  to  a simple  matrix  equation  which 
includes  unknown  charges  onlv  over  each  conductor.  The  number  of  unknowns 
is  considerably  reduced  over  that  of  method  //l,  but  the  variety  of  problems 
which  can  be  treated  is  more  limited.  The  second  method  also  may  be  used 
to  obtain  the  values  of  discontinuity  capacitances  in  microstrip  [56] , [65] . 

Fig.  4-14  shows  computed  values  of  discontinuity  capacitance  for  a junction 
between  two  microstrip  lines  of  different  strip  widths. 

The  third  method  utilizes  an  integral  solution  of  the  microstrip  problem 
and  is  applicable  to  multiple  dielectric  and  multiple  ground  planes  [70],  [76]. 
Computation  time  is  relatively  lengthy.  Fig.  4-15  shows  the  capacitance  of 
a three- layer  microstrip  computed  by  this  third  method.  More  complete  informa- 
tion on  microstrip  applications  is  covered  in  references  [40],  [56],  [651,  [70], 
[76],  and  in  volume  #6  of  the  RADC  Method  of  Moments  Applications  Series  [146]. 


4.5  APERTURE  COUPLING 


Electromagnetic  coupling  through  an  aperture  in  a plane  screen  can  be 

* 


treated  by  the  method  of  moments.  The  aperture  is  replaced,  using  Babinet's 
principle,  with  the  complementary  metal  disk,  which  is  in  turn  replaced  by 
a wire  grid  of  appropriate  size  and  shape.  The  fields  in  the  aperture,  and 
near  and  far  fields  beyond  the  aperture  may  then  be  computed  by  the  method 


* 

The  method  of  moments  as  used  here  allows  one  to  treat  apertures  up  to  about 
one  square  wavelength  in  area. 
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Figure  4-14a  A sudden  change  in  width  of  microstrip. 
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of  moments  [3.47].  Consider  a A/2  by  A/4  rectangular  aperture  in  a ground 
plane.  This  aperture  is  replaced  by  a complementary  disk  which  is  wire- 
gridded  as  shown  in  Fig.  4-16.  A plane  wave  polarized  in  the  v-direction 
is  incident  normally  upon  the  aperture.  Incident  field  intensity  is  one 
volt/meter  and  frequency  is  100  MHz.  Fig.  4-16  shows  resultant  disk  cr-rent 
I as  computed  by  two  different  method-of-moments  techniques.  Fig.  4-17 
shows  a comparison  of  computed  and  experimental  data  for  the  y-component  of 
the  resultant  aperture  field.  Note  that  the  agreement  is  very  good  near  the 
central  portion  of  the  aperture.  Near  the  edges  the  agreement  is,  of  course, 
less  good  since  fields  approach  zero  or  infinity  (E  ® at  y ° 0.125A  and 
+ 0 at  x = 0.25A)  and  measurement  problems  arise.  Fig.  4-18  shows  a 
near  field  plot  of  the  electric  field  in  the  shadow  region  beyond  the  aper- 
ture along  a circle  of  radius  0.5A  lying  in  the  x-z  plane  (see  Fig.  4-16). 

A plane  wave  polarized  in  the  y-direction  is  normally  incident  upon  the 
aperture.  The  experimental  data  shown  in  Figs.  4-16  to  4-18  was  taken  at 
Sandia  Corporation. 

In  the  method-of-moments  technique  for  aperture  coupling,  arbitrary 
angles  of  incidence  and  arbitrary  polarization  can  be  readily  treated. 
Examples  are  given  in  reference  [147],  In  addition,  coupling  to  wires 
beyond  the  aperture  may  also  be  treated.  More  complete  information  on 
aperture  coupling  by  the  method  of  moments  is  available  in  references  [147], 
[148]. 
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Figure  4-16  Equivalent  disk  currents  on  a rectangular  grid,  due  to  a plane  wave 
normally  incident  upon  a rectangular  aperture  in  a plane  screen, 
polarized  parallel  to  the  short  side.  (E1  * 1 volt/meter,  f » 100  mhz). 
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Figure  4-17  Fields  in  a rectangular  aperture  in  a plane  screen,  due  to  a 

plane  wave  normally  ircident,  polarized  parallel  to  the  short  side. 
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V - CONCLUSIONS 

The  method  of  moments  has  proven  to  be  a successful  method  for  treating 
a wide  variety  of  problems,  many  of  which  could  not  be  treated  heretofore 
with  a high  degree  of  accuracy.  Examples  of  such  problems  are  coupling  in 
the  presence  of  parasitics,  near  field  and  radiation  hazard  prediction, 
transmission  line  interference,  etc.  The  basis  of  its  application  to 
thin-wire  structures  is  described  in  this  volume.  The  method  is  extremely 
flexible  and  conceptually  simple,  and  it  is  applicable  in  theory  to  most 
problems  ir.  elect  ;omagnetic  theory.  Illustrati/e  examples  of  its  appli- 
cation are  presented. 
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APPENDIX  A - JUSTIFICATION  OF  THE  7HIN-WIRE  MODEL 


In  this  appendix,  Che  effects  of  finite  radius  are  considered  and  an 

attac.pt  is  made  to  justify  the  filamencary  model  used  in  section  1.2.  Fig. 

A-l  shows  a fat  dipole  with  surface  currents  J , J, , J . There  are  several 

z 9 P 

steps  that  need  to  be  justified:  (1)  neglect  of  circumferential  currents, 

J.,  (2)  neglect  of  radial  currents,  J , (3)  neglect  of  4>  variations  in 
<p  P 

longitudinal  currents  J , and  ;4)  replacement  of  z-directed  surface  currents 

Z 

with  a filamentary  model. 

The  thin-wire  model  may  be  explained  in  several  different  ways.  The 
results  of  some  classical  scattering  problem);  are  useful.  Consider  a plane 
wave  incident  upon  a conducting  cylinder  of  radius  "a"  and  Infinite  length 
(Figure  A-2)  with  its  axis  parallel  to  the  z-direction.  The  incident  electric 
field  has  intensity  Eq  and  is  polarized  in  either  the  z direction  (parallel 
to  the  cylinder  axis)  or  in  the  y direction  (perpendicular  to  the  cylinder 
axis),  as  shown  in  Fig.  A-2.  The  solution  of  the  resultant  boundary  value 
problems  for  the  currents  and  .ittered  fields  are  given  below  [149): 


ELECTRIC  FIELD  PARALLEL  TO  THE  CYLINDER  AXIS* 


Surface  Currants 


rn  Jn* 


_r  H<2)(ka) 

r.»-r  n 


(A-l) 


t (J) 

In  this  appendix,  J (x.  (x)  represent  Bessel  functions  of  the  first 

kind  and, Henkel  functions  n of  the  second  kind,  respectively,  and  J*(x) 
and  H*^'(x)  rtprasant  their  derivatives  with  respect  to  x.  n represents 
the  order  of  the  respective  functions,  y Is  Euler's  constant  (1.741). 
k la  tha  fraa  spacs  wava  number.  (p,4,t>  are  the  cylindrical  coori lnates. 
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(2) 

A study  of  Che  small  argument  forms  of  H (ka)  shows  that  the  n *■  0 term 

n 

of  the  Infinite  series  dominates  for  small  ka  and  thus 
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% 


H(2)(ka) 

o 


1 - ~ log(^p) 


-1 2 


log (.8905  ka) 


where  log  indicates  the  natural  logarithm. 

JE 


j . - — 2. I 

z toga  log (.8905  ka) 


(ka  <<  1) 


-12tt  Eq 

I (total  current)  = 2?raJ  n : n-v'r  ", — r 

z z uu  log (.8905  ka) 


for  (ka  <<  1) , 


(A-2) 


(ka  <<  1)  (A-3) 


The  n = 0 term  of  the  current,  which  becomes  the  total  current  for  small  ka, 
varies  slowly  [l/(log(. 8905  ka))]  with  ka  for  small  ka. 

Scattered  Fields 


Near  Fields 
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2kp  log (.8905  ka) 


(ka  <<  1) 


(A-7) 


Thus,  as  ka  -*■  0,  the  ^-independent  currents  dominate  and  the  radiated  fields 


become  ^-independent.  Note  especially  that  the  near  electric  field  is 
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proportional  to  H^  (kp) /log (.8905  ka)  (see  Eq.  (A-5)).  The  term  H^(fcp) 
o o 

is  proportional  to  the  tield  of  a filament  along  the  z-axis,  whereas  the 
term  l/log(.8905  ka)  corresponds  to  the  variation  of  total  longitudinal 
surface  current  for  small  ka  (See  Eq.(A-3)).  Thus,  for  small  ka,  the  ra- 
diated fields  of  a filament  are  equal  to  those  of  the  distributed  surface 
current  model  with  the  same  total  current,  for  p > a.  This  is  the  basis 
for  step  (4),  the  replacement  of  the  surface  model  with  a filamentary  model. 


II.  ELECTRIC  FIELD  PERPENDICULAR  TO  THE  CYLINDER  AXIS 


Surface  Currents 


<p  wka  L H( 
n**-«  n 


J-n  ejn* 

„(2)’ 


(A-8) 


where  Hq  ■ E^/q  “ incident  magnetic  intensity. 

The  current  induced  has  a dipole  (cos  <£)  term  at  wel?  as  a ^-independent 
(n  ® 0)  term.  Both  are  important  for  small  ka.  The  <>- independent  current 
represents  electrically  small  loops  of  electric  current  (for  small  ka) . Elec- 
trically small  loops  are  equivalent  to  magnetic  dipoles  and  thus  the  <J>- 
independent  electric  current  is  equivalent  to  a magnetic  current  filament 
of  infinite  length. 


Scattered  Fields 


Far  Fields 


2 | « J^(ka)  e“ 

"l  „(2)?(ks 
n°-»  n 


(A-9) 
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as  ka  ♦ 0 


Note  that  both  the  dipole  (cos  4>)  and  the  magnetic  filament  contribute  as 

v 2 

ka  -*■  0.  However,  the  scattered  fields  are  proportional  to  (ka)  in  contrast 
to  the  scattered  fields  for  the  parallel  case  which  are  proportional  to 
1/ (log. 89  ka.)  . In  other  words,  a thin  wire  may  have  both  longitudinal  and 
circumferential  currents.  As  ka  becomes  small,  the  fields  radiated  by  the 
circumferential  component  decrease  rapidlv  with  ka,  whereas  the  fields  radi- 
ated by  the  longitudinal  component  decrease  slowly  with  ka.  This  corresponds 
to  the  well-known  fact  that  the  back-scattering  cross  section  of  a thin  wire 
parallel  to  the  incident  field  varies  slowly  with  wire  radius  and  that  it  is 
much  greater  than  that  of  a thin-wire  perpendicular  to  the  incident  field. 

We  can  combine  the  results  of  the  two  problems  to  express  the  scattering 
due  to  arbitrary  incident  fields  with  electric  field  components  both  parallel 
and  perpendicular  to  the  cylinder  axis. 

Note  that,  for  arbitrary  polarization  of  the  incident  fields,  both 
longitudinal  and  circumferential  currents  are  excited.  The  fields  radiated 
by  the  longitudinal  currents  will  always  dominate  as  ka  ^ 0,  except  for  the 
special  case  where  there  is  no  component  of  parallel  to  the  cylinder. 

Furthermore,  note  that  because  of  the  dominance  of  the  n «*  0 term  in  (A-l) 
and  (A-4) , the  variations  in  <t>  of  the  longitudinal  current  are  unimportant. 

Thus,  for  plane  waves  normally  incident  upon  infinite  cylinders  of  small 
ka,  the  scattered  fields  radiated  by  the  longitudinal  currents  dominate  over 
thope  radiated  by  the  circumferential  currents.  Furthermore,  the  fields 
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radiated  by  the  ^-independent,  longitudinal  currents  dominate  over  those  ra- 
diated by  the  ^-dependent  longitudinal  currents.  The  above  statements  are 
valid  in  both  the  near  and  far  fields.  The  onlv  exception  to  the  above 
statements  is  the  degenerate  case  where  the  incident  field  is  exactly  per- 
pendicular to  the  cylinder  axis. 

In  order  to  apply  the  above  classical  analysis  to  antenna  problems  in- 
volving thin  wires,  one  must  generalize  the  analysis  in  several  respects.  One 
must  consider  non-nonnal  (oblique)  incidence,  near  field  sources  rather  than 
plane  waves,  wires  of  finite  rather  than  infinite  length.  Suffice  it  to  say 
that  none  of  these  generalizations  completely  invalidate  the  conclusions 
reached  above.  However,  some  of  the  conclusions  are  weakened  somewhat.  For 
instance,  consider  two  parallel  wires  extremely  close  together  (a  fraction  of 
a wire  radius).  One  may  be  considered  as  a near  field  source  for  the  other. 

It  is  possible,  for  instance,  that  a current-carrying  wire  extremely  close 
to  a cylinder  may  excite  significant  ^-dependent  longitudinal  currents.  Even 
though  these  currents  radiate  less  effectively  than  the  ^-independent  longi- 
tudinal currents,  there  are  undoubtedly  cases  where  the  radiated  fields  of 
the  ^-dependent  currents  are  significant.  For  instance.  If  two  wires,  both 
electrically  small  in  radius,  are  placed  in  extremely  close  proximity  to  each 
other,  then  one  would  expect  that  (a)  the  currents  would  concentrate  near  the 
portions  of  the  wires  closest  to  each  other  and  (b)  the  <J>— de;  endent  longitu- 
dinal current  would  become  significant.  How  close  must  the  wires  be  to  in- 
validate the  assumptions  used  in  typical  thin-wire  programs?  The  results  of 
a static  analyris  provide  scr..e  guidelines.  Fig.  A-3  shows  two  parallel  cylin- 
ders of  infinite  length.  A voltage  is  applied  between  the  cylinders.  The 

79 


I 


kaa.. . . . srfasa^a  ijfaB  ttktasaaAae 


charge  distribution  on  the  cylinders  is  desired.  This  problem  may,  of 
course,  be  solved  by  classical  methods.  For  more  than  two  cylinders,  or 
for  dielectric-coated  cylinders,  the  method  of  moments  can  be  used. 

Clements  (150)  has  treated  such  problems  by  the  method  of  moments  using 
sinusoidal  expansion  functions  (i.e.,  a Fourier  series  decomposition  in  <}>) 
for  the  induced  charge  distribution.  The  induced  charge  is  represented  as 
follows : 

N 

j(4>)  ■ 2 cos  n(£  + s*n  n$ 

n=0 

For  the  case  shown  in  Fig.  A-3,  the  origin  of  the  coordinate  system 
may  be  chosen  such  that  all  the  B^'s  are  zero.  Fig.  A-4  shows  the  ratio 
Ai/Aq,  i.e.,  the  ratio  of  the  coefficient  of  the  dipole  term  to  that  of  the 
(^-independent  term,  as  a function  of  D/2a.  Figs.  A-5a,b,c  show  the  relative 
magnitudes  of  the  Fourier  coefficients  for  the  cases  D/2a  * 1.05,  1.5,  5.0 
respectively.  Note  that  for  very  close  spacing  such  as  D/2a  = 1.05  (.1  radius 
separation),  the  dipole  coefficient  exceeds  that  of  the  ^-independent  term. 

For  D/2a  =1.5  (one  radius  separation  between  cylinders)  the  dipole  coefficient 
is  comparable  to  that  of  the  (^-independent  term.  For  D/2a  = 5.0  (4  diameters 
separation),  the  magnitude  of  the  dipole  coefficient  is  less  than  1/6  that 
of  the  ^-independent  term. 

The  TEM  transmission  line  problem  is,  of  course,  directly  related  to  the 

static  problem,  and,  as  a result,  the  ratios  A^Aq  of  Figs.  A-5a,b,c  are 

interpreted  directly  as  the  Fourier  coefficients  of  the  longitudinal  surface 

current  density  Jz  of  the  transmission  line  problem  where: 

N 

J * 7 A cos  n<}>'  + B sin  n$ 

z L n n 
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Figure  A-5a  Fourier  coefficients  of  the  two-wire  problem  (D/2a  = 1.05) 
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From  this  analysis  then,  one  would  suspect  that  the  ^-variation  of  the  longi- 
tudinal current  could  not  be  neglected  unless  the  wires  were  several  diameters 
apart. 

In  addition  to  the  ^-directed  and  <p—  dependent  currents  one  must  justify 
the  neglect  of  the  p-directed  current.  Consider  the  four  disks  of  p-directed 
currents  (two  at  the  gap  and  two  at  the  ends  of  the  dipole  of  Fig.  A-l) . For 
small  ka,  the  radial  (p-directed)  currents  are  ineffective  radiators  for  several 
reasons.  All  of  the  current  elements  are  electrically  short  and  are  ineffective 
radiators  for  that  reason.  Furthermore,  the  directions  of  the  currents  on 
opposite  sides  of  a disk  are  such  as  lead  to  cancellation.  Further,  the 
currents  on  tho  two  gap  disks  tend  to  cancel.  The  currents  on  the  end  disks  are 
considerably  smaller  than  those  on  the  gap  disks.  Note  that  very  close  proximity 
of  parallel  wires  would  upset  some  but  not  all  of  these  conclusions.  Thus, 
in  most  problems,  there  is  a high  order  of  cancellation  of  the  effects  of 
the  radial  currents  and  they  may,  therefore,  be  neglected. 

In  summary,  for  thin-wire  antennas  (ka  <<  1),  the  surface  currents  J. 

<P 

and  J may  be  neglected.  The  circumferential  variation  of  J may  be  neglected 

P 2 

and  the  resultant  uniform  surface  current  J may  be  replaced  by  a filament  of 
current  along  the  axis  of  the  wire.  These  conclusions  are  based  upon  a considera- 
tion of  plane  wave  scattering  by  a cylinder.  For  thin  wires  very  close  to- 
gether (within  a few  radii  of  each  other),  these  assumptions  may  be  invalid. 
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APPENDIX  B - Z^  FOR  STRAIGHT  PARALLEL  WIRES  WITH  PULSE 
EXPANSION  FUNCTIONS  AND  IMPULSIVE  WEIGHTS 


In  this  appendix  z^  is  expressed  In  terms  of  the  ij>  function.  Fig. 

B-la  she  .•»  the  typical  subsections  A£.,  Li  .'ith  currents  i . , i..  z,, 

,v  i 3 i j ij 

may  be  expressed  as  follows  (see  Eq.  (1-26)). 

-iMr^r’l 

7-1 TT  d2'  (B-l) 

I L,  ~ L I 


'ij 


V ’ J 

i Li 


j 


A£iA£j  ,a2 


* + k >* 


where  tf>  *•  -rr—  --  i TT  dz* 

A£  J ( .r I r.  - r.  j 

J AO  A 


(B-2) 


Note  that  if/  is  proportional  to  the  tangential  magnetic  vector  potential  A 

z 

at  one  subsection  (AJ^)  due  to  unit  current  on  another  (A£^).  Fig.  B-la 
shows  the  basic  geometry. 

The  origin  may  be  translated  to  the  center  of  subsection  Li ^ as  in 
Fig.  B-lb  without  loss  of  generality.  For  simplification,  coordinates 
^xi’^i,zP  C^e  point  i are  replaced  in  this  appendix  with  (x,y,z), 

(r,0,<f>),  (p»$»z)  in  rectangular,  spherical,  and  cylindrical  coordinates  as 
shown  in  Fig.  B-lb.  Because  of  rotational  symmetry,  if>  is  Independent  of  $ 
and  depends  only  on  p,z,  and  the  electrical  length  (kAfc^)  of  subsection  A£^ 
Eq.  (B-2)  becomes 


ij 


jawAfcjAS,^  - 


Li^Ll 

jue 


iii 


3z 


(B-3) 


upon  replacing  zJ  with  z. 
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Figure  B-l  Typical  parallel  subsections  for  thin-wire  antenna  problems 
(a)  general  location  (b)  translation  to  the  origin. 


^^s*i^;38eg!^5^5!w«a>5*5ga^!5^^g585RS0&«9^;^s,,,»s^^®w 


3<^/3z  mav  be  approximated  by  finite  differences  as  follows: 

3(J)  <p(x,v,z  + ^|)  - i|/(x,y,2  - “) 

3z  " Az 


(B-4) 


where  Az  is  a small  length  increment  - not  necessarily  equal  to  A£^  and 
where  4i(x,y,z  + ~)  is  proportional  to  the  potential  at  point  (x,y»z  + ^) 
with  the  source  subsection  A £j  fixed  as  in  Fig.  (B-lb)  with  center  at  the 
origin.  Equation  (B-4)  is  written  in  shorthand  notation  as 


3*  *<*  ft  ~ *<-  *§> 

3z  Az 


(B-5) 
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Note  that  3i^/3z  may  also  be  evaluated  by  keeping  the  field  point  fixed  and 


shifting  the  source  subsection  Ai^ , ie 

3,  -rv(+^)  - i'-*£n 

3z  Az' 


(B-6) 


where  ¥(+  Az')  is  proportional  to  the  potential  at  field  point  i with  the  center 


a y 1 

of  the  source  subsection  (A£^)  shifted  to  (0,0,  — ) i.e.,  shifted  upwards 


by  an  increment 


A z * Az 1 

~Y~  • V ( 2~  ) indicates  a downward  shift. 


Eqs.  (B-5)  and  (B-6)  are  related  to  the  identity 


3r 

3z 


3r 

3z’ 


where  r 
Now 


J(x-x')2  + (z-z’)*  + (v-y*) 
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(B— 7) 
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S Az 1 

- ['?( r~)'J  represents  the  derivative  of  ^ with  a shifted  source  sub- 

OZ  2 

section;  this  may  be  evaluated  by  shifting  the  field  point  by  Az/'2: 

. . n-  ^ , + 4|>  - , - % 

fi  [*«=!*)  i 5 — Hr* — - — - 


(B— 8) 


h !',‘±r)1 


Az1  Az 


f(+  ~ + ^f)  - ?(+  ~ ) 


(B-9) 


where  ¥(+  ”2“* represents  $ with  A£^  of  Fig.  B-lb  shifted  — j—  upwards 


and  A2.J  shifted  — downwards,  etc. 


Collecting  terras: 

-A"  - - TT7TT  W+  H1.  + 


Az.  , Az*  Az. 

— ) - H+  — - — ) 


1/  Az*  . Az.  (/  Az*  Az.  . 
- i>(~  ~ » + “j)  + <K 2“. j)] 


(B-10) 


Substituting  in  Hq.  (B-3)  yields 


Af,.A^ 


Az*  . Az.  , .Az*  Az. 


• j— vy + eSp1  t*<+  *r-  + 2> 

- *1- *£,+*$**(■  *£, 


(B-ll) 


Note  that  Eq.  (B-ll)  involves  the  evaluation  of  five  different  functions 

with  vectors  r,.,  r,.(+  +) , r, .(+  -),  etc.  as  shown  in  Fig.  B-2,  and  corre- 
—ij  “ij  ~ ij 

spending  terras  0^(4-  ")»  2i j vhere  z^C^-  “)  represents  the  distance 

from  shifted  field  point  i to  the  perpendicular  bisector  plane  of  shifted 
source  subsection  Ai^,  etc.  Each  of  the  different  <p  functions  may  be 
evaluated  by  the  methods  outlined  by  Harrington  [1],  [2]. 
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Eq.  (15-11)  has  been  derived  under  a consideration  of  current  distribu- 
tion alone;  the  charge  distribution  is  determined  by  the  equation  of  continu- 
ity. Thus  pulse  current  expansion  functions  lead  to  impulsive  charges  at 
the  ends  of  each  subsection. 

it  is  interesting  to  compare  the  foregoing  result  based  on  Eq.  (1-19) 
and  assuming  pulse  currents  and  impulsive  charges,  with  one  based  on  Eq.(l-12) 
and  derived  under  the  assumption  of  pulse  expansion  functions  for  both 
currents  and  charge  [3].  This  latter  combination  is  inconsistent  mathematically 
but  represents  a more  realistic  charge  distribution. 

If  Az  » A£^,  Az‘  « A£ j , then  Eq.  (B-ll)  becomes 


$ 

>5 


z^  «*  jwuA^Afc.^  + [t|<(+  +)  - &(+  -)  - <p(-  +)  + <p( )]  (B-12) 

which  ag.rees  with  Eq.  (22)  of  (3].  Thus  the  two  methods  (pulse  currents 

and  impulsive  charges  - pulse  currents  and  pulse  charges)  lead  to  identical 

results  if  the  increments  Az,  Az*  used  for  finite  differences  are  half 

2 2. 

subsections. 
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APPENDIX  C - SYMMETRY  ABOUT  THE  MAIN  DIAGONAL  FOK  THE  MATRICES  [?.)  AND  [D] 


In  general,  for  the  method  of  moments,  the  matrices  [z]  and  [D]  are 
not  symmetric  about  the  main  diagonal,  i.e. 


2u  * zji 


Du  * V 


v/here  z^  is  the  voltage  at  suosecticn  1 due  to  unit  current  density  on 


subsection  j (in  time-harmonic  problems)  and  D^.  is  the  potential  at  sub- 


section i due  to  unit  charge/unit  charge  density  on  subsection  j [40]. 

If  the  weighting  and  expansion  functions  are  identical  (Galerkin's  method, 
a variational  solution),  then  the  matrices  [z]  or  [D]  are  symmetrical 
about  the  main  diagonal. 

For  parallel  current  segments  A£.,,  A2,.  of  equai  length,  z 


i’  “~j  ~ “ij  ji 

For  parallel  segments  of  unequal  length,  or  for  skewed  segments,  is  not 


necessarily  equal  to  z^.  Similar  distinctions  apply  to  charge  filaments 


A2^,  Afcj  and  functions. 


Consider  an  example  from  electrostatics  consisting  of  two  perpendicu- 
lar, thin,  charged  wires,  two  typical  subsections  cf  which  are  hown  in 
Fig.  C-l.  Subsections  AS^  and  A£j  are  shown  with  line  charge  densities  X i 


and  Xj  and  subsection  centers  i and  j respectively.  Subsection  lengths  differ: 


A i - c 


and  the  separation  of  subsections  is  equal  to  "a* 


is  the  potential  at 


D , may  be  defined  as  potential  due  to  unit  charge  or  unit  charge  density. 
The  former  definition  is  used  here. 
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point  i due  to  unit  charge  on  subsection  ASU  and  D ^ is  the  potential  at 
point  j due  to  unit  charge  on  subsection  A£j  . 


a+b  . , , 

1 f XJ  dy 


Dij  * (Lb)  “ ^ V 4 ire  J y'  “ 4ire  b lo8(1  + a} 
J 1 1 o 1 } o 


(C-l) 


where  is  the  potential  at  point  i due  to  A , 


£ c_  . , , K 2 

_ ax'  1 2 J 4 + (a  2)  (C-2) 

Jji  4ttc  c j / ““  4ttc  c °8  ~ l"- 2 ~ 7 9 

J (x  , + (a  + j'  " 2'  + > — + (a  + 2> 


where  log  indicates  the  natural  logarithm. 
If  b = c,  then  Eq.  (C-2)  becomes 


Dji  4neob  *°8 


| 4-  J,"»ab+f 

- \ + Ja2  + ab  + ■— 


(C-3)  | 


Clearly  is  not  equal  to  D^,  even  for  b = c.  However,  as  the  separa- 
tion ncr eases  (b/a  and  c/a  approach  zero),  approaches  D^: 


’ll  ‘ 4^b  108  (1  + b/a)  4^7  t“/*  ~ » 

J r>  /> 


(C-4) 


4*e  c 
o 


c + /T,  , b2+c2  , b' 

1 , , * 4a £ X 

= + /[“'7777m 


I’1 

-!  + a<1  + fe> 


1 , 1 ■<•  c/(2a  t b) 

47tc  c 108  x - c/ (2a  + b> 
o 


^~log  [1  + 2c/ (2a  + b)  j 


1 

4rtt  a 


(b/a  <<  1,  c/a  <<  1) 
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Fig.  C-2  shows  a plot  of  R (where  R ■ as  a function  of  b/a, 

for  the.  special  case  b a c.  Note  that  the  matrix  [D]  Is  highly  asymmetric 
about  the  main  diagonal  for  large  b/a,  and  that  1). ^ is  indistinguishable 
from  for  b/a  < 0.3.  Stewart  [128]  has  shown  explicitly  that,  if 
pulse  weighting  and  expansion  functions  are  used  in  this  problem,  " Dji‘ 
This  is  true  in  general;  the  Galerkin  method  always  yields  matrices  sym- 
metric about  the  main  diagonal. 
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